Mapping the three-body system - decay time and reversibility

In this paper we carry out a quantitative analysis of the three-body systems and map them as a function of decaying time and intial conguration, look at this problem as an example of a simple deterministic system, and ask to what extent the orbits are really predictable. We have investigated the behavior of about 200 000 general Newtonian three body systems using the simplest initial conditions. Within our resolution these cover all the possible states where the objects are initially at rest and have no angular momentum. We have determined the decay time-scales of the triple systems and show that the distribution of this parameter is fractal in appearance. Some areas that appear stable on large scales exhibit very narrow strips of instability and the overall pattern, dominated by resonances, reminds us of a traditional Maasai warrior shield. Also an attempt is made to recover the original starting conguration of the three bodies by backward integration. We find there are instances where the evolution to the future and to the past lead to different orbits, in spite of time symmetric initial conditions. This implies that even in simple deterministic systems there exists an Arrow of Time.Comment: 8 pages, 9 figures. Accepted for publication in MNRAS. Includes low-resolution figures. High-resolution figures are available as PNG

A Hybrid N-Body Code Incorporating Algorithmic Regularization and Post-Newtonian Forces

We describe a novel N-body code designed for simulations of the central regions of galaxies containing massive black holes. The code incorporates Mikkola's 'algorithmic' chain regularization scheme including post-Newtonian terms up to PN2.5 order. Stars moving beyond the chain are advanced using a fourth-order integrator with forces computed on a GRAPE board. Performance tests confirm that the hybrid code achieves better energy conservation, in less elapsed time, than the standard scheme and that it reproduces the orbits of stars tightly bound to the black hole with high precision. The hybrid code is applied to two sample problems: the effect of finite-N gravitational fluctuations on the orbits of the S-stars; and inspiral of an intermediate-mass black hole into the galactic center.Comment: 12 pages, 15 figures, accepted for publication in MNRA

Asteroids in the Inner Solar System I - Existence

Ensembles of in-plane and inclined orbits in the vicinity of the Lagrange points of the terrestrial planets are integrated for up to 100 million years. The integrations incorporate the gravitational effects of Sun and the eight planets (Pluto is neglected). Mercury is the least likely planet, as it is unable to retain tadpole orbits over 100 million year timescales. Both Venus and the Earth are much more promising, as they possess rich families of stable tadpole and horseshoe orbits. Our survey of Trojans in the orbital plane of Venus is undertaken for 25 million years. Some 40% of the survivors are on tadpole orbits. For the Earth, the integrations are pursued for 50 million years. The stable zones in the orbital plane are larger for the Earth than for Venus, but fewer of the survivors are tadpoles. Both Venus and the Earth also have regions in which inclined test particles can endure near the Lagrange points. For Venus, only test particles close to the orbital plane are stable. For the Earth, there are two bands of stability, one at low inclinations (i < 16 degrees) and one at moderate inclinations (between 24 degrees and 34 degrees). The inclined test particles that evade close encounters are primarily moving on tadpole orbits. Our survey of in-plane test particles near the Martian Lagrange points shows no survivors after 60 million years. Low inclination test particles do not persist, as their inclinations are quickly increased until the effects of a secular resonance with Jupiter cause de-stabilisation. Numerical integrations of inclined test particles for timespans of 25 million years show stable zones for inclinations between 14 and 40 degrees.Comment: 20 pages, 21 figures, Monthly Notices (in press

Chaos in the one-dimensional gravitational three-body problem

We have investigated the appearance of chaos in the 1-dimensional Newtonian gravitational three-body system (three masses on a line with $-1/r$ pairwise potential). We have concentrated in particular on how the behavior changes when the relative masses of the three bodies change (with negative total energy). For two mass choices we have calculated 18000 full orbits (with initial states on a $100\times 180$ lattice on the Poincar\'e section) and obtained dwell time distributions. For 105 mass choices we have calculated Poincar\'e maps for $10\times 18$ starting points. Our results show that the Poincar\'e section (and hence the phase space) divides into three well defined regions with orbits of different characteristics: 1) There is a region of fast scattering, with a minimum of pairwise collisions and smooth dependence on initial values. 2) In the chaotic scattering region the interaction times are longer, and both the interaction time and the final state depend sensitively on the starting point on the Poincar\'e section. For both 1) and 2) the initial and final states consists of a binary + single particle. 3) The third region consists of quasiperiodic orbits where the three masses are bound together forever. At the center of the quasiperiodic region there is the periodic Schubart orbit, whose stability turns out to correlate strongly with the global behavior.Comment: 24 pages of text (REVTEX 3.0) + 21 pages of figures. Figures are only available in paper form, ask for a preprint from the author

Collisional dynamics around binary black holes in galactic centers

We follow the sinking of two massive black holes in a spherical stellar system where the black holes become bound under the influence of dynamical friction. Once bound, the binary hardens by three-body encounters with surrounding stars. We find that the binary wanders inside the core, providing an enhanced supply of reaction partners for the hardening. The binary evolves into a highly eccentric orbit leading to coalescence well beyond a Hubble time. These are the first results from a hybrid ``self consistent field'' (SCF) and direct Aarseth N-body integrator (NBODY6), which combines the advantages of the direct force calculation with the efficiency of the field method. The code is designed for use on parallel architectures and is therefore applicable to collisional N-body integrations with extraordinarily large particle numbers (> 10^5). This creates the possibility of simulating the dynamics of both globular clusters with realistic collisional relaxation and stellar systems surrounding supermassive black holes in galactic nuclei.Comment: 38 pages, 13 figures, submitted to ApJ, accepted, revised text and added figure

Long-Term Evolution of Massive Black Hole Binaries. III. Binary Evolution in Collisional Nuclei

[Abridged] In galactic nuclei with sufficiently short relaxation times, binary supermassive black holes can evolve beyond their stalling radii via continued interaction with stars. We study this "collisional" evolutionary regime using both fully self-consistent N-body integrations and approximate Fokker-Planck models. The N-body integrations employ particle numbers up to 0.26M and a direct-summation potential solver; close interactions involving the binary are treated using a new implementation of the Mikkola-Aarseth chain regularization algorithm. Even at these large values of N, two-body scattering occurs at high enough rates in the simulations that they can not be simply scaled to the large-N regime of real galaxies. The Fokker-Planck model is used to bridge this gap; it includes, for the first time, binary-induced changes in the stellar density and potential. The Fokker-Planck model is shown to accurately reproduce the results of the N-body integrations, and is then extended to the much larger N regime of real galaxies. Analytic expressions are derived that accurately reproduce the time dependence of the binary semi-major axis as predicted by the Fokker-Planck model. Gravitational wave coalescence is shown to occur in <10 Gyr in nuclei with velocity dispersions below about 80 km/s. Formation of a core results from a competition between ejection of stars by the binary and re-supply of depleted orbits via two-body scattering. Mass deficits as large as ~4 times the binary mass are produced before coalescence. After the two black holes coalesce, a Bahcall-Wolf cusp appears around the single hole in one relaxation time, resulting in a nuclear density profile consisting of a flat core with an inner, compact cluster, similar to what is observed at the centers of low-luminosity spheroids.Comment: 21 page

Dynamics of two planets in co-orbital motion

We study the stability regions and families of periodic orbits of two planets locked in a co-orbital configuration. We consider different ratios of planetary masses and orbital eccentricities, also we assume that both planets share the same orbital plane. Initially we perform numerical simulations over a grid of osculating initial conditions to map the regions of stable/chaotic motion and identify equilibrium solutions. These results are later analyzed in more detail using a semi-analytical model. Apart from the well known quasi-satellite (QS) orbits and the classical equilibrium Lagrangian points L4 and L5, we also find a new regime of asymmetric periodic solutions. For low eccentricities these are located at $(\sigma,\Delta\omega) = (\pm 60\deg, \mp 120\deg)$, where \sigma is the difference in mean longitudes and \Delta\omega is the difference in longitudes of pericenter. The position of these Anti-Lagrangian solutions changes with the mass ratio and the orbital eccentricities, and are found for eccentricities as high as ~ 0.7. Finally, we also applied a slow mass variation to one of the planets, and analyzed its effect on an initially asymmetric periodic orbit. We found that the resonant solution is preserved as long as the mass variation is adiabatic, with practically no change in the equilibrium values of the angles.Comment: 9 pages, 11 figure

Interaction of massive black hole binaries with their stellar environment: II. Loss-cone depletion and binary orbital decay

We study the long-term evolution of massive black hole binaries (MBHBs) at the centers of galaxies using detailed scattering experiments to solve the full three-body problem. Ambient stars drawn from a isotropic Maxwellian distribution unbound to the binary are ejected by the gravitational slingshot. We construct a minimal, hybrid model for the depletion of the loss cone and the orbital decay of the binary, and show that secondary slingshots - stars returning on small impact parameter orbits to have a second super-elastic scattering with the MBHB - may considerably help the shrinking of the pair in the case of large binary mass ratios. In the absence of loss-cone refilling by two-body relaxation or other processes, the mass ejected before the stalling of a MBHB is half the binary reduced mass. About 50% of the ejected stars are expelled ejected in a "burst" lasting ~1E4 yrs M_6^1/4, where M_6 is the binary mass in units of 1E6 Msun. The loss cone is completely emptied in a few bulge crossing timescales, 1E7 yrs M_6^1/4. Even in the absence of two-body relaxation or gas dynamical processes, unequal mass and/or eccentric binaries with M_6 >0.1 can shrink to the gravitational wave emission regime in less than a Hubble time, and are therefore "safe" targets for the planned Laser Interferometer Space Antenna (LISA).Comment: Minor revision. 10 pages, 7 figures, ApJ in pres